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ACCOMPLISHMENTS: 


A series of impact experiments on anhydrite CaSO-tin which vaporized sample 
accelerates an element in a velocity interferometer, generate velocity data that we have 
recently reanalyzed using an explicit entropy generating finite difference code. The shock 
pressure required from the onset, and complete vaporization of 30% porous (19 ± 2.8 and 
98.7 ±11.3 GPa) and 70% crystal density anhydrite is 52 ± 3 and 122 ±13 GPa. (These 
initial data were later revised as explained below). Using observed acid leaching in non- 
marine K/T ejecta in N. America, and the sharp global increase in 87 Sr/ 86 Sr ratios 
recorded at 65 Ma in marine rocks, we demonstrated that global acidification is primarily 
due to the SO2 released by anhydrite volatilization, and not HNO3 formed from bolide- 
induced air pyrolysis. Shock temperatures for crystal CaCC>3 are measured from 3000 to 
7000 K in the 90 to 160 GPa pressure range. These temperatures are much lower than 
calculated theoretically indicating that possibly bond breakdown at the shock front is 
occurring. This is the first mineral in which this effect has ever been seen. New data 
defining the ion species which are produced upon impact of volatilization of metals and 
minerals using a pulse UV laser to simulate intense shock heating from a projectile 
impact indicate that in shock experiments we can for the first time study the speciation of 
neutrals using a moderate resolution (0.2 dalton) time-of-flight mass spectrometer. 
Measurements of the gas species from a series of proposed impact experiments appear to 
be quite feasible. We will attempt these experiments in the next year. Measurements of 
the impact induced shock wave decay in SiC>2 and GeC>2 glass are underway to measure 
these pure oxide properties. Predictive calculations indicate that the pressure will decay 
as r 27 in the phase transition regime, versus a much lower rate of r' 1 8 , if a phase 
transition does not occur. 

We have successfully developed the techniques necessary to conduct impact 
experiments on ice at very low temperatures. We employ the method of embedding 
gauges within a target to measure the shock wave and material properties. This means 
that our data are not model dependent; we directly measure the essential parameters 
needed for numerical simulations of impact cratering. Our preliminary experiments are 
published in Stewart and Ahrens [2000]. Since then we have developed a new method 
for temperature control of icy targets that ensures temperature equilibrium throughout a 
porous target. Graduate student, Sarah Stewart-Mukhopadhyay, is leading the work on 
ices and porous materials as the main thrust of her thesis research. 

The work on ice and ice-silicate mixtures follows our previous work: “Low- 
velocity impact craters in ice and ice-saturated sand with implications for martian crater 
count ages” Croft et al. [1979], “Identification of ice VI on the Hugoniot of ice Ih” 
Gaffney and Ahrens [1980], “Fragmentation of ice by low velocity impact” Lange and 
Ahrens [1981], “The dynamic tensile strength of ice and ice-silicate mixtures” Lange and 
Ahrens [1983], and “Impact experiments in low-temperature ice” Lange and Ahrens 
[1987]. Our previous work has focused on icy materials with no porosity, and we 
propose to extend our research to include porous ice and porous ice-silicate mixtures. 
There is little shockwave data for porous ice [Bakanova et al., 1976; Furnish and 
Boslough, 1996], and none of the data was acquired under conditions applicable to the 
outer solar system. The solid ice Hugoniot is only defined for initial temperatures above 
-20C [Gaffney, 1985], 
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Our program uniquely measures the properties of ice at temperatures directly 
applicable to the solar system. Previous experiments were conducted at ambient 
temperatures soon after removing the target from a cold environment, usually just below 
freezing, or in a room just below freezing (e.g. Arakawa [1999]). Since ice has an 
extremely complicated phase diagram, see Figure 1, it is important to conduct 
experiments at lower temperatures to determine the true outcome of impacts in the outer 
solar system (e.g. Gaffney and Matson [1980]). 

Our work has focused on the inherent material properties by measuring the shock 
wave directly; this complements the macroscopic observations and immediately provides 
the parameters necessary to extend this research to the gravity regime. 

Our numerical simulations of impacts in porous ice under very low gravity 
conditions, such as found on comets, show that the final crater size and shape is very 
dependent on the dynamic strength of the material. There is some controversy over 
whether or not solid ice has a Hugoniot Elastic Limit (HEL), the maximum dynamic 
stress the material may bear before yielding. The data, based primarily on experiments 
near -10C, have alternatively been interpreted as a phase change (melting) instead of an 
HEL. This demonstrates the importance of conducting experiments at the appropriate 
temperature to understand the ice phases involved. Until our research, there has never 
been a measurement of the HEL in porous ice and it was assumed that porous ice had no 
inherent dynamic strength. We have discovered a small HEL in porous ice. We are 
conducting further experiments to verify this detection and understand its porosity and 
temperature dependence as well as the effects of thermal processing (sintering). 
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(b) (after Larson [1984]) 

FIGURE 1. Importance of initial temperature on 
shock Hugoniot of H 2 Q. (a) The different segments of 


3 






Pressure (GPa) 

(c) H 2 0 Phase Diagram [ Lobban et al. , 1998] 


the ice Hugoniot centered at - IOC correspond to the 
shocked phase as illustrated in (b). (b) An estimate of 
the path of the ice Hugoniot shown in (a), (c) H 2 0 
phase diagram down to absolute zero. The 
importance of the initial temperature is evident from 
the many different phases of ice that can be reached 
by moderate shock pressures. 


We have conducted more recovery experiments to measure sound velocity deficits 
in shock-damaged materials, and we are working on relating the deficits to the impact 
conditions: impact velocity and projectile mass. 

Ultrasonic velocities are a good indicator of impact damage because the wave 
speed is a strong function of the crack or pore density in the material (e.g. O’Connell and 
Budiansky [1977]). 

Figure 2. P-wave velocities in San Marcos gabbro. 
The sound speed is measured in targets impacted 
with a spherical Al projectile at the noted velocity. 
The data shown are from directly beneath the center 
of the crater. The sound speed deficit near the crater 
indicates impact damage. The higher the impact 
velocity, the lower the wave speed, and the larger the 
impact damage. The regions further from the crater 
remain undamaged and retain the pre-impact wave 
speeds. The difference in undamaged sound 
velocities between the samples is attributed to 
heterogeneities in San Marcos gabbro. 3a error 
bars. 
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There is currently very little work in this area in the geophysics community 
whereas headway has been made in the ceramics community (e.g. Camacho and Ortiz 
[1996]). Our previous work includes: “Mechanical properties of shock-damaged rocks” 

He and Ahrens [1994] and “Stress wave attenuation in shock damaged rock” Liu and 
Ahrens [1997]. Our examination of the literature on terrestrial impact craters (e.g. 

Pilkington and Grieve [1992]) shows that there seems to be a log-log relationship 
between the crater diameter and the depth of the impact-damaged zone [Ahrens et al., 

2001]. We are quantifying this relationship based on the impact physics and the new 
insight gained from the research on impact damage in ceramics. Our program on impact 
damage in geologic materials will allow the planetary science community to quantify the 
effects of cratering in another dimension: below the surface. We also expect to conduct 
multiple impact experiments to understand cratering in previously damaged materials. 



4 





Figure 3. Damaged zones beneath terrestrial 
craters inferred from seismic refraction 
measurements. Using published data, we have 
estimated the relationship between the crater 
diameter and the depth of the damaged zone 
for a dozen terrestrial craters. This empirical 
relationship has not yet been explained 
quantitatively. The least squares fit is 
log H (km) = -0.367 + 0.673 log D (km) where 
D is the crater diameter and H is the depth of 
the damaged zone. 
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We have revisited the question of how much vapor was produced from the K/T 
impact [Gupta et al, 2001]. New analysis of data on shock vaporization of anhydrite 
[Yang and Ahrens, 1998] has produced revised values for the shock pressures required for 
incipient and complete vaporization. Together with a revised size for the transient 
diameter of the Chicxulub crater [O'Keefe and Ahrens, 1999], the amount of sulfuric 
vapor that could have produced sulfuric acid aerosols has been decreased to 1 1-39 
gigatons, compared to 200 gigatons estimated by Pope et al. [1997], Assuming a 
maximum production of sulfuric acid aerosols, the global temperature decreased 12-19K 
for about a decade, compared to the 5-31 K drop for 12 years reported by Pope et al. 
[1997]. Since the Earth’s surface temperature during the K/T epoch was 18-20K warmer 
than present values (e.g. Savin and Yeh [1981]), this cooling event could not produce 
global freezing conditions at the Earth’s surface. 

Even with this new result, the global effects of the K/T impact cannot be 
accurately assessed without an understanding of the speciation in the vapor plume. If 
most of the sulfur ends up in SO 2 , the global impact may still be considerable and long- 
lived. We have been developing a method to detect impact vapor species using a time-of- 
flight mass spectrometer (TOFMS). 

Another application of impact volatilization is the accretion of the Earth and the 
role of serpentine, which is illustrated by the following reactions: 

Mg 3 Si 2 05(0H)4 — > Mg 2 Si04 + MgSiOj + 2H 2 0 
(serpentine) (forsterite) (enstatite) (water) 

( 1 ) 

MgSi0 3 + Fe + H 2 0 -> ~ Fe 2 Si0 4 Mg 2 Si0 4 +H 2 

z z* 

(enstatite)(iron)(water) (fayalite) (forsterite) (hydrogen) 

( 2 ) 

Under this grant, we conducted a detailed study of reaction (1), which included 
determining the Hugoniot of serpentine (e.g. Tyburczy [2001]; Tyburczy et al. [1991]; 


5 




Tyburczy [1990], Ahrens [1990]). We also note that serpentine is one of the most 
abundant hydrous minerals in primitive meteorites [Lange et al., 1985; Ringwood, 1979]. 
It has been inferred that serpentine was abundant in planetesimals and was the carrier of 
most of the water to the terrestrial planets. Abe and Matsui [1985] directly utilized the 
Caltech impact devolatilization data on serpentine to define their model of a super 
greenhouse overlying a magma ocean during accretion of the Earth, a condition that 
prevailed for a large fraction of the accumulation of the planet. After about 0.26 of the 
final mass of the Earth had accreted, further impacts devolatilized creating a thick 
hydrosphere. 
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